Ab initio study of pressure stabilised NiTi allotropes: pressure-induced 
transformations and hysteresis loops 
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Changes in stoichiometric NiTi allotropes induced by hydrostatic pressure have been studied 
employing density functional theory. By modelling the pressure-induced transitions in a way that 
imitates quasi-static pressure changes, we show that the experimentally observed B19' phase is (in 
its bulk form) unstable with respect to another monoclinic phase, B19". The lower symmetry of 
the B19" phase leads to unique atomic trajectories of Ti and Ni atoms (that do not share a single 
crystallographic plane) during the pressure-induced phase transition. This uniqueness of atomic 
trajectories is considered a necessary condition for the shape memory ability. The forward and 
reverse pressure-induced transition B19'-f->-B19" exhibits a hysteresis that is shown to originate from 
hitherto unexpected complexity of the Born-Oppenheimer energy surface. 

PACS numbers: 61.50.Ah, 61.50.Ks, 62.20.fg, 64.30.Ef, 64.60.My, 81.30.Hd 
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I. INTRODUCTION 

Nickel-titanium alloys belong to the important class 
of shape-memory materials^— . Their properties include 
super-elasticity, excellent mechanical strength and ductil- 
ity, good corrosion resistance and bio-compatibility (im- 
portant for example in medical applications), and high 
specific electric resistance (allowing the material to be 
easily heated by an electric current). The shape memory 
effect is governed by a martensitic transformation from 
a high-temperature austenitic phase (cubic B2, CsCl- 
structure) into a low-temperature martensitic phase. X- 
ray experiments on single crystals^ and neutron mea- 
surements on powder samples^ revealed the low temper- 
ature phase to be a monoclinic B19' structure (see Fig.[TJ 
7 « 97.8°) with P2i/m space group. In addition, a rhom- 
bohcdral R-phase^ with P3 space group was found dur- 
ing multi-step martensitic transformations^— under the 
following conditions: (i) off-stoichiometric composition, 
(ii) presence of substitutional or interstitial impurities, 
and/or (iii) formation of precipitate phases. 

Several theoretical studies on the low temperature 
martensitic phase of stoichiometric NiTi alloys have been 
performed. The intense search has been motivated in 
part by the fact that theoretically predicted structures 
do not unambiguously agree with those detected exper- 
imentally. For example Huang, Ackland, and Rabe— 
concluded that the B19' structure is unstable with re- 
spect to a higher-symmetry base-centered orthorhombic 
(BCO, in some studies also termed B33) structure (see 
Fig. [TJ 7 w 107°). These conclusions were based on 
systematically cross-checking several distinct DFT meth- 
ods, functionals, and implementations (FLAPW, PAW, 
USPP, GGA, LDA, ABINIT, VASP, etc.). The analysis 
also considered a carefully selected shear transformation 
path connecting all three structures B2 (7 = 90°), B19', 



and BCO, since they are characterized by a specific value 
of the crystallographic angle 7. Very similar results were 
reported by Wagner and Windl— and by Guda Vishnu 
and Strachan— . The latter authors^ also predicted a 
new phase (B19") characterized by 7 m 102.5° and with 
practically identical energy to the BCO phase. Finally, a 
barrier-less transformation path between the B2 and the 
BCO phases as a sequence of several special deformation 
modes was demonstrated in Ref. 

[13 

Various explanations of the discrepancy between (i) the 
apparent stability of the B19' phase as observed in low- 
temperature experiments and (ii) the instability of the 




FIG. 1. Atomic geometry of the investigated B19'-like phases. 
The various structures considered in this study alternate in 
the lattice parameters a, b, c, monoclinic angle 7, and in- 
ternal positions (see text for details). Larger blue spheres 
correspond to Ni, smaller gray spheres to Ti atoms. The high- 
lighted planes are used to characterize the structural ability to 
accommodate the shape-memory effect (see Sec. IIIIDj) . The 
picture was generated using the VESTA package^. 
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B19' phase predicted by theoretical calculations (for 
T = K) have been proposed: Recent theoretical works 
of Sestak, Cerny, and Pokluda^ and Zhong, Gall, and 
Zhu— suggest that the B19' may be stabilized by the 
presence of (nano)twins that are often experimentally 
observed^. As another possibility Huang, Ackland, and 
Rabe— suggested that the B19' structure could be sta- 
bilized by residual stresses that are frequently present in 
experimental samples. Since the equilibrium volume is 
predicted to be smaller for the B19' structure than for 
the BCO phased, one may expect the BCO structure to 
transform into the B19' phase under compressive loads. 

Considering this variety of mechanisms active in NiTi 
and in order to understand how external strains effect 
the stability of the various phases, we systematically ex- 
plore the potential energy surface (PES). To comple- 
ment previous studies, we focus solely on martensitic 
phase transformations induced by volumetric changes, 
i.e., hydrostatic pressure. Our choice is motivated by 
the fact that (i) stress/strain fields in NiTi alter process- 
parameters of the martensitic transformations (such as 
e.g., the transition temperature) and (ii) these actual 
stresses and strains in experimental samples are difficult 
to measure and are often not known. Focusing on volu- 
metric changes, we show an unexpectedly complex PES. 
This complexity results in transformation mechanisms 
that exhibit hysteresis effects not reported in previous 
studies. From a methodological point of view, we also 
show that it is difficult to include internal variables ex- 
plicitly in the PES since they are responsible for mctasta- 
bility of and the newly discovered hysteresis processes. 



II. COMPUTATIONAL DETAILS 

The calculations were performed using density func- 
tional theory (DFT) 2 ^^ in the generalized gradient ap- 
proximation (GGA-PBE'96)^ 2 . as implemented in the Vi- 
enna Ab-initio Simulation Package (VASP)22i2i. All 
monoclinic structures were studied using four-atom cells 
with different external and internal parameters, while a 
two-atom cell was used for the B2 phase. As the total en- 
ergy differences among different phases are rather small, 
it was necessary to ensure convergence of the energy be- 
low f meV per formula unit (f.u.), i.e., one Ni and one Ti 
atom. Therefore, the plane wave cutoff energy was set to 
400 cV and a 24 x 16 x 18 fe-point Monkhorst-Pack mesh 
was used to sample the Brillouin zone of the monoclinic 
allotropes studied. 



A. Computational Methodology: Quasi-Static 
Volumetric Changes 

The computational approach usually employed for 
studying the effect of hydrostatic pressure is based on 
determining the total energy as function of volume. The 
hydrostatic pressure in the system is obtained by fitting 



the equation of stated to the calculated energy-volume 
data points. Because the B19' and the BCO phases are 
structurally similar and differ only slightly in few internal 
(atomic coordinates) and external (lattice constants and 
the angle 7) parameters, the multi-dimensional Born- 
Oppenhcimcr potential energy surface (PES) is expected 
to be quite complex, exhibiting many local minima. In 
order to explore the impact of hydrostatic pressure on 
phase stability and martensitic phase transformations 
among different NiTi allotropes, we determined the PES 
as function of (i) the atomic volume, (ii) Ni atom x- 
axis internal coordinate, and monoclinic angle 7 (see de- 
tails below). In order to systematically map the com- 
plex PES, we adopted a quasi-static (QS) approach, 
within which the volume is increased/decreased gradu- 
ally in an adiabatic-like maimer (see detailed explana- 
tion in Appendix |A"[) . This not usually used approach 
allows for more realistic simulations of gradually increas- 
ing/decreasing pressures since it closely imitates experi- 
mental conditions. 

III. RESULTS AND DISCUSSION 

A. The monoclinic allotropes under hydrostatic 
load 

The QS simulations were initiated using the previ- 
ously identified ground states for each phase (B19' and 
BCO). Subsequently, both structures were to evolve 
quasi-statically under applied volumetric changes. Fig. [2] 
summarizes results from four separate simulations. From 
the initial configuration (either B19' or BCO), we first 
gradually increased the volume to the maximum value 
studied here (blue circles and dashed lines in Fig. [2]), and 
subsequently decreased the volume in the QS manner to 
the lowest calculated value (blue circles and solid lines 
in Fig. [2]). Similarly, we proceeded in the opposite direc- 
tion: from the initial state we first decrease the volume 
down to the minimum investigated value (red triangles 
and dashed lines in Fig. [2]), and then increased it to the 
maximum (red triangles and solid lines in Fig.[2|). We ap- 
plied these two forward-and-backward runs to both the 
B19' and BCO starting configurations. 

When starting the QS volumetric changes with the 
BCO phase, the angle 7 ranges between 105° for 
20A 3 /f-u. (« 60GPa, compression) and 108° for 
33 A /f.u. (« -20GPa, expansion). The internal coor- 
dinate XNi remains almost constant at the value w 0.915. 
In contrast to what was suggested by Huang, Ackland, 
and Rabe—, no transition to the B19' phase is observed 
within this fairly broad range of hydrostatic pressures. 

A very different behavior is obtained for the B19' start- 
ing configuration. The application of positive hydrostatic 
pressures (red dashed path in Fig. [2]) first changes the 
starting angle 7 abruptly from « 100° to rj 94°. Further 
decreasing the volume results in only small changes of 
the angle 7. Again, no transition to the BCO structure 
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FIG. 2. Calculated dependence of (a) the monoclinic angle, 7, and (b) the internal coordinate, xm, on the volume obtained 
from the QS simulations (see text). Along the B19' compression/decompression paths a region when an additional phase (B19") 
becomes stable is found. Three ranges of monoclinic angle 7, that correspond (from the lowest values of 7 to the highest) to 
B19', B19" and BCO phases, respectively, are separated by horizontal dashed lines. 



is predicted. Surprisingly, when negative pressures are 
applied (volumetric increase, see the dashed blue path in 
Fig- E]), the angle 7 changes to approximately 103°. The 
resulting unit cell geometry and the internal coordinates 
no longer correspond to values typical for either the B19' 
or the BCO state. A similar behavior is demonstrated 
in Fig. [2Jd for the volumetric dependence of the internal 
coordinate xm- The structural parameters of this state 
are very similar to the B19" phase described by Guda 
Vishnu and Strachan— . 

Our results allow to disregard early suggestions of a 
B19'-O-BC0 transition induced by hydrostatic pressure. 
Rather, we conclude that hydrostatic pressure, similar 
to shear deformations^, transforms B19' into B19". In 
contrast to monoclinic shear—, hydrostatic strain does 
not drive a transition towards BCO. Finally, we find that 
the BCO phase is stable with respect to the hydrostatic 
deformations and does not transform to B19' (or B19"). 

B. Origin of the B19'oB19" hysteresis 

A closer look at the reaction pathways in Fig. [3] reveals 
the presence of a narrow hysteresis loop. To explain its 



origin, we have analyzed the PES along both transition 
paths (resulting from increasing and decreasing volume) . 
We expressed the total energy as a function of a single 
external parameter, volume V, and one selected internal 
parameter, here the x Nl position. We chose the latter 
parameter because, unlike the angle 7, it can easily be 
kept constant in available DFT implementation and it 
provides clear ranges defining the two phases, B19' and 
B19" (see Fig. Hand Table HJ). 

Using these two parameters we have calculated the 
total potential energy surface £? PES (x Nl , V) (Fig. (3^)- 
As expected, the B19" structure is a stable phase (a 

global minimum at V « 27.5 A /f.u. and a?Ni ~ 0.94), 
while B19' is not associated with any minimum, indi- 
cating that in a fully relaxed environment (hydrostatic 
pressure p = 0) this phase is unstable. To investigate 
the influence of external strain we consider the enthalpy 
H(x m , V lP ) = E PES (x m , V) + pV with p being the (hy- 
drostatic) external applied pressure. Increasing the pres- 
sure p shifts the equilibrium volume of the B19" phase to- 
wards smaller values (Fig.[3p). In addition, at sufficiently 
high pressures, a new minimum occurs that represents 
for p = 5.22 GPa the B19' phase. Fig. [3}d explains also 
neatly the occurrence of the hysteresis. To go from one 
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FIG. 3. Potential enthalpy surface for p = for fully relaxed states (a) and for p — 5.22 GPa (b) along the reaction paths volume 
V and Ni internal atomic position x Kl . The red, blue and violet symbols mark different volumetric loading conditions and 
corresponding trajectories (indicated by arrows). Part (c) shows vertical const ant- volume cuts of the p — 5.22 GPa enthalpy. 



phase to the other, even at the critical pressure (p = 5.22 
GPa) where both phases have identical enthalpy, a bar- 
rier along the constant-volume paths exists. Since in an 
adiabatic transformation only the nearest local minimum 
is reachable, the trajectory follows the original path even 
though this minimum is no longer the energetically most 
favorable one. 

To demonstrate this further, we plot in Fig. (3fc the 
energy profiles at fixed volumes (vertical profiles corre- 
sponding to the PES in Fig.[3Ji). The figure clearly shows 
that upon increasing the volume (i.e., following the red 
triangles), the structure is trapped in a local energy val- 
ley, and transforms to B19" only when the energy barrier 
completely flattens. A similar mechanism happens also 



in the opposite direction (i.e., following the blue circles). 

To further confirm this hypothesis, we performed an 
additional test. We started from a B19"-like structure, 
but from a volume (~ 27.3 A /f.u.) only slightly larger 
than that at which the B19'— >B19" transition occurs 
(ra 26.9 A /f.u.). This pathway is marked by purple stars 
in Fig. [3l This pathway also crosses the B19" minimum 
and eventually joins the B19"— s>B19' branch of the origi- 
nal hysteresis, i.e., the one corresponding to volume com- 
pression ("blue circle" data points). 

From Fig. [3)d we can further deduce that on chang- 
ing the transformation coordinates from volume to x Nl 
(which is closely related to the monoclinic angle 7) qual- 
itatively different paths result. In this scenario only a 
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single minimum for a fixed value of x the transformation 
coordinate (the horizontal cuts of the PES) is obtained, 
instead of the two-minima for the vertical cuts shown in 
Fig. Consequently, in this case corresponding to a 
shear mode transformation along the angle 7, no hys- 
teresis occurs. 

We thus conclude that the structural complexity of 
the B19' and B19" phases and (related to it) the 
multi-minimum character of the PES are the origin for 
the transformation hysteresis under hydrostatic loading. 
This is in contrast to the structurally much more distinct 
phases, B2 and B19, as shown by Kibey, Sehitoglu, and 
Johnson—. 

In contrast to what may be expected the volume in- 
creasing (red triangles) and decreasing (blue circles) data 
points in Fig. [3] do not coincide in the region away from 
the hysteresis loop. The reason is the energy difference 
between the states (expanding (red triangle) and shrink- 
ing (blue circle) volume) at a constant volume is in the 
order of (or smaller than) lmeV/f.u.. This value is be- 
low the numerical accuracy of the present calculations. 
The apparent discrepancy is thus simply a consequence 
of extremely flat valleys of the PES corresponding to the 
B19' and, in particular, B19" phases. Increasing the cal- 
culation accuracy (albeit at significant CPU costs) is ex- 
pected to result in a closer correspondence of the two 
pathways. 

C. Structural parameters 

The potential energy surface shown in Fig. [3^ provides 
sets of quasi-static energy-volume data. These data sets 
that can be individually analyzed using the Murnaghan 
equation of stated. Following this approach we get the 
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true ground state properties of all phases and can assign 
a pressure value to each data point. Part of these results 
are presented in Fig. HJ In this graph all data points 
are plotted, i.e., from both branches of the hysteresis in 
Fig. [2j As a criterion for separating the B19' and B19" 
phases we used the internal coordinate a; Nl : a structure 
with x Nl > 0.97 is B19'-like otherwise it corresponds to 
the B19" phase (sec Fig. [2b). 

Focusing on the most interesting region close to the 
equilibrium volumes (Fig. 2]), we could have easily mis- 
taken the BCO and BI9" states as a single phase if we 
had not performed a thorough analysis of internal coordi- 
nates and lattice parameters. The bulk moduli and their 
pressure derivatives, as well as the equilibrium volumes 
and the structure energy differences from the Murnaghan 
equations of state are summarized in Table [H together 
with all the equilibrium structural parameters. As can 
be seen, the energy of the BI9" state is equal to that of 
the previously predicted BCO state within the numerical 
accuracy of our calculations. 

A comparison of the structural parameters of our BI9" 
phase with those obtained by Guda Vishnu and Stra- 
chan— reveals some differences, the largest being in the 
monoclinic angle 7 (100.0° predicted here vs. 102.4° cal- 
culated by Guda Vishnu and Strachan— ). These are 
likely to be consequences of using different deformation 
modes (hydrostatic vs. shear). Despite these small dif- 
ferences we regard these two structures as "flavors of the 
same phase and thus use the same name B19" for both 
of them. 

Finally, differences between our structural parameters 
and those reported in the literatu re 15 ' 17 of the B19' phase 
stem from the fact that in the earlier studies the mono- 
clinic angle 7 was fixed to the experimental value (« 98°) 
while we allowed for a full structural relaxation. As men- 
tioned in the previous section, performing a full relax- 
ation reveals that at ambient pressure the B19' phase is 
unstable with respect to the B19" phase. 

We note that because some phases are stable only in 
certain volume (pressure) range, the E(V) data points 
could not be computed over the whole volumetric range. 
For example, the B19' phase data points are only avail- 
able for volumes smaller than w 27 A /f.u.. The pre- 
dicted properties of the ground states should nevertheless 
be reasonably accurate as the number of data points ob- 
tained is sufficient to perform a numerically robust fitting 
to the equation of state. 

An advantage of the QS approach is that the energy- 
volume data points are less scattered (i.e., less influenced 
by the complexity of the NiTi PES) and their numerical 
analysis is therefore more robust. Consequently, the ini- 
tial states (minima of the non-QS energy-volume curves) 
differ from final states (minima of the quasi-static E(V) 
curves). An example is shown in Fig. [5] and illustrates 
these differences: the non-QS value for the monoclinic 
angle of the B19' structure is 7 100.5° while the QS 
analysis gives 7 w 94.5°. Since we consider the QS cal- 
culations (that mimic the experimental pressure increase 
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TABLE I. Calculated external and internal parameters of the B2, B19', BCO, and B19" structures in their ground states. 
The atomic volumes V cq , bulk moduli Bo, their pressure derivatives B' , the total energy differences with respect to the BCO 
phase of the equilibrium states AE = E C q — E^° , lattice parameters a, b, c, and the angle 7 together with the internal x 
and y coordinates of the Ni and Ti atoms predicted in the present study, and compared with the literature data of Huang, 
Ackland, and Rabe— and Guda Vishnu and Strachan 17 . All the extensive quantities, such as the total energy differences and 
equilibrium volumes, are listed per two-atom NiTi formula unit (f.u.). 



or decrease) more accurately, the discrepancy between 
the QS and non-QS ground states demonstrates the ne- 
cessity to compute the energy-volume dependence quasi- 
statically. 

In order to determine the critical pressure needed for 
the B19'-(->-B19" phase transition, we calculated the en- 
thalpies H of both phases (Fig. [5^). Since the enthalpy- 
pressure data calculated for different phases are similar, 
an analytical formula for the enthalpy functional was 
used for the fitting. 

Fig. [SJd shows the enthalpy difference between a given 
phase and the phase with the lower enthalpy for a specific 
pressure. The corresponding transition pressure from 
B19" to B19' is 5.22 GPa. We note that this pressure 
dramatically reduces when applying the analytical for- 
mula by a linear (0.47 GPa) or quadratic (2.01 GPa) fit. 
This finding clearly demonstrates the necessity of using 
the analytical expression from Ref. [2?] based on the Mur- 
naghan equation of stated. It should be further noted 
that this single- value transition pressure neglects kinetic 
effects. 

The value of the critical hydrostatic pressure, ~ 5 GPa, 
above which B19' becomes stable may be compared with 
the value of 1 GPa when applying shear stress^. Ex- 
ploring the complexity of possible mechanisms active in 
NiTi alloys, our study and that by Wagner and Windl— 
also complement recent work by Sestak, Cerny, and Pok- 
ludaJ^ proposing a twinning mechanism for the stabiliza- 
tion. 



D. Ability of the monoclinic allotropes to show a 
shape memory effect 

In contrast to the orthorhombic BCO phase (which 
has a too high symmetry to account for the shape 
memory effect^ ), both the B19' and B19" structures 
possess only a lower (monoclinic) symmetry. The 
lower symmetry guarantees that the atomic austenitc- 
martensitc transition pathway within the unit cell is 
unique. Thus, the structural phase can, in principle, 
store the shape information since all the atoms remain 
situated in the Erickscn-Pittcri neighborhood of their 
austcnitc counterparts^. 

In order to recover the BCO lattice from the B19" 
structure, the Ni atom above the base center has to move 
into the plane defined by two Ni atoms from the unit cell 
basal plane and one of the neighboring Ti atoms (see 
Fig. [!}. The internal atomic positions then fulfill a spe- 
cific geometric relation that can be quantified by a struc- 
tural parameter S. Employing the internal structural pa- 
rameters fri, ^Ni, (Ti, and ^Ni as defined in Fig. [TJ the 
parameter 8 is given as: 

6=p/p±. (1) 

When 5=1, the corresponding Ni and Ti atoms are 
located within a single plane and the structure is too 
symmetric to keep an atom-to-atom relationship with the 
B2 lattice necessary for the shape-memory effect. As 
the S parameter distinguishes whether these Ni and Ti 
atoms are, or are not, in a planar arrangement, it will be 
termed "planarity" parameter: A deviation from 5 = 1 
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FIG. 5. (a) The theoretically predicted enthalpy of the B19' 
and B19" phases over the whole range of studied pressures, 
(b) The differences with respect to the phase that minimizes 
the enthalpy for a given pressure. 



is a prerequisite to store the shape information. 

The planarity parameter 5 for both the BCO and B19" 
phases is shown in Fig. [SJ Due to the high symmetry en- 
vironment, the Ti atoms in the BCO phase (full circles 
in Fig. [6]) are very stable in their location in spite of 
high hydrostatic pressures (small volumes). 5 is a con- 
stant function of volume for the BCO phase indicating 
that the high symmetry is preserved during the hydro- 
static loading. In contrast, the volume dependence of 
the planarity parameter of B19" (empty circles in Fig. [5]) 
deviates from 5—1 over the whole range of volumes 




volume (A 3 /f.u.) 



FIG. 6. The planarity parameter, 5, describing the internal 
symmetry of the phases as a function of volume. The value 
of 1 found for the BCO phase indicates that the Ti atoms 
are located within the same atomic plane as the Ni atoms 
(see Fig. []} and the internal geometry of the phase is too 
symmetrical. The B19" values deviating from 5=1 indicate 
an ability to store the shape information. 



studied. This is in agreement with the analysis of the 
B19" structure by Guda Vishnu and Strachan— . Within 
the theory of symmetry-dictated extremal - — , the BCO 
phase represents a structure with a symmetry-dictated 
energy minimum. 



IV. CONCLUSIONS 

We report on first-principles calculations of pressure- 
induced transitions in stoichiometric NiTi allotropes. 
Complementing previous studies that focused on shear 
strains and twinning mechanisms, we have systemati- 
cally explored the complex potential energy surface of 
NiTi under well-defined generic volumetric changes. We 
kept the volume constant at each simulation step of the 
pressure-induced transitions and relaxed all other struc- 
tural degrees of freedom with respect to the total energy. 
By repeating these steps in a quasi-static manner, we 
closely mimicked experimental conditions. 

In contrast to previous theoretical studies of shear 
deformations^—, the BCO phase does not transform 
into the experimentally observed B19' phase when apply- 
ing hydrostatic pressures. We ascribe this stability of the 
BCO allotrope to the high symmetry of this structure. In 
contrast, the B19' structure distorts under pressure into 
another, newly identified, monoclinic phase, B19". This 
phase is located structurally in between the B19' and 
BCO phases. We find that the B19" phase has an energy 
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comparable to that of the BCO phase. 

The complexity of the Born-Oppcnheimer potential en- 
ergy surface results in a pressure-induced B19'<->B19" 
transition that exhibits a previously unreported hystere- 
sis. The latter could be related to the inherent multi- 
dimensional nature of the potential energy surface in 
NiTi. The B19" structure has a lower symmetry than 
the BCO phase. As a consequence, the B19" structure 
can be the basis of the shape memory effect. 
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Appendix A: The quasi-static approach 

In the following we explain the essential difference be- 
tween quasi-static (QS) and non-QS calculations. The 
more common non-QS approach can be regarded as a se- 
ries of energy-volume data points obtained by (i) start- 
ing with an identical internal atomic coordinates and the 
overall cell shape and (ii) changing only the overall vol- 
ume of the unit cell. A so-called relaxed state, i.e., a set 
of external, {£"}, and internal parameters, {C™}, for a 
given volume, V n , is then found by minimizing the total 
energy, E, as a function of both internal and external 
parameters (except for the volume): 

non-QS: = mk E(V = V n , {&, CD 

where & = £°,Ci = Cf reflects the fact that the start- 
ing configuration for all E(V) data points (labeled with 
n) is the same. The relaxed states found in the non-QS 
approach can exhibit phases different from the starting 
one if pressure- induced structural transitions occur in the 



studied system. The non-QS computational approach 
may fail to properly simulate experimental conditions as 
the non-QS states that were obtained by discontinuously 
changed volume may differ from those found in experi- 
ments (in which volume and/or strain are always varied 
continuously). The advantage of the non-QS simulation 
is that all the calculations can be performed indepen- 
dently in a parallel manner, i.e., calculations for different 
volumes can be distributed over all available computa- 
tional units (processors). 

In contrast to the non-quasi-static simulations, the 
quasi-static (QS) simulations can not be parallelized as 
they proceed by a subsequent series of calculations con- 
sisting of the following steps. First, all parameters are 
energy-relaxed for a certain starting state that frequently 
corresponds to the equilibrium conditions, i.e., zero hy- 
drostatic pressure. Then, a small change of the volume is 
applied and new sets of internal and external parameters 
are obtained by the total energy minimization with re- 
spect to the structural parameters (volume being fixed). 
With these new relaxed parameters (i) a small volumet- 
ric change and (ii) subsequent structural optimization are 
repeated. The two steps are then repeated so as to cover 
the whole range of volumes: 

QS: min E(V = V n ,£i, &) -> {# , CD ■ 

£»=£»"- \C*=C 

The QS simulation mimics a compression of the struc- 
ture in case of negative volumetric changes and decom- 
pression in case of positive ones. The denser is the 
mesh of the calculated volumes the better correspon- 
dence should be achieved with the experimental com- 
pression/ decompression processes. The QS procedure en- 
sures that the system evolves smoothly from one local 
minimum into another and the pressure-induced transi- 
tions including the phase transition path in a complex 
configurational space can be studied. In contrast, a non- 
QS search algorithm may results in (non-physically) dis- 
continuous jumps in the atomic trajectories. 
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